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We analyze how to exploit Brillouin scattering for the purpose of cooling opto-mechanical devices 
and present a quantum- mechanical theory for Brillouin cooling. Our analysis shows that significant 
cooling ratios can be obtained with standard experimental parameters. A further improvement of 
cooling efficiency is possible by increasing the dissipation of the optical anti-Stokes resonance. 



Introduction. — A lesser known quality of Brillouin 
scattering is its ability to scatter photons in the anti- 
Stokes direction [IJ, while cooling the corresponding me- 
chanical vibration. The recent experimental observation 
of such Brillouin cooling [2] suggests that a model is re- 
quired to describe the potential of this new system for 
ground-state cooling [3j [4] . In such light-sound interac- 
tions, light is scattered to both Stokes and anti-Stokes 
side bands, which heat and cool the system (Fig. ^p). 

In order to break the material heating-cooling sym- 
metry, we study here a setup involving a resonator with 
an asymmetric resonance structure. Figure [T^i exempli- 
fies our proposal to resonantly enhance the anti-Stokes 
process for cooling the mode, while at the same time 
off-resonantly attenuating the Stokes process to prevent 
heating. Obtaining a resonator that is proper for Bril- 
louin cooling is challenging. This is because two opti- 
cal resonances that have almost the same optical fre- 
quency but different propagation constants are needed 
in order to conseerve both the energy and momentum 
which are given to light by the acoustical phonon. One 
type of cavity that allows such optical-resonance pairs is 
whispering-gallery mode resonators [5j EJ. In such res- 
onators, the transverse (radial-polar) order of one mode 
can compensate for the frequency difference originating 
from the non-similar longitudinal (azimuthal) order of 
the other. This provides a pair of modes with differ- 
ent azimuthal wavevectors, but nearby frequencies, as 
experimentally observed via the resulting stationary in- 
terference pattern [5j [6]. The energy flow in Brillouin 
anti-Stokes cooling [2 that is analyzed here, is opposite 
in respect to the Stokes excitation process [7H9]. Addi- 
tionally, cooling here is a spontaneous process. 

Derivation of Optoacoustic Coupling Rate. — First we 
derive the optoacoustic coupling rate T opt associated with 
Brillouin cooling. Consider a purely longitudinal density 
fluctuation that is propagating along L through a cross- 
sectional area A. A displacement field u causes an elastic 
potential of \AT(d z u) 2 , where T is the spring constant 
per area. Additionally, the kinetic energy will be \Apix 2 
where p is the mass density. The interaction with the 
light field E is given by \A^[d z u)E 2 , where 7 = p^ 
is the electrostrictive constant which relates a change in 
density to a change in permittivity. Thus our Hamilto- 
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FIG. 1. a) Sample resonator spectrum showing resonant en- 
hancement possible at both pump and anti-Stokes frequencies, 
with attenuation at Stokes frequency, b) Energy diagram for 
Brillouin heating and cooling processes. 
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where the omitted parts refer to the Hamiltonian of the 
electromagnetic field inside the medium. Here tt(z) = 
Apu is the momentum density. 

Quantization of the Electromagnetic Field and Sound 
Wave. — For the electromagnetic field, consider a single 
polarization subjected to periodic boundaries. The elec- 
tric field is of the form: 



E(z) = J2 E k [a k e ikz + h.c] 



(2) 



where E k is the zero-point amplitude of the electric field. 
Knowing the total energy of the electric and magnetic 
fields in free space, E 2 and B 2 , to be per mode for 
the ground-state, we solve for the zero point fluctuation 
of the electric field, 



E k = y/(hw k )/(2ALe) . 



(3) 



Similar to the electric field, the longitudinal sound wave 
will be quantized in the following form: 



H z ) = ^2 u k \b k e lkz + h.c 



(4) 



Substituting Eq. Q into Eq. ([!]) and solving for the 
coefficient of b k b\,, we obtain the expressions: 



u k = y / (h)/(2ALpuj k ) , 



(5) 
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where uo k = \k\v s is the acoustic dispersion relation, and 
v s = y/T/p is the speed of sound. 

Acousto- optical interaction. — Having derived expres- 
sions for both E k and u k we turn our attention to the 
optomechanical interaction term. Substituting expres- 
sions for E k and Uk from Eqs. Q and Q into and 
using the rotating wave approximation, the Hamiltonian 
simplifies to the following form: 
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= h J29k, q (K + bl q )(al+ q a k + a k a{ +q ) , (6) 

k,q 

where = (1/2) AL^qu q E k E k+q . 

The pump mode (p) is coupled to another optical mode 
(aS) and an acoustical mode in which energy conservation 
(u p = uj a s — ft) and momentum conservation (k p = k a s — 
q) are fulfilled [9 . Here, q and ft represent the wavevector 
and frequency of the acoustical mode. Conservation of 
both quantities is possible in the kind of experimental 
setup studied here due to high order transverse optical 
modes [5j [6|. Taking into account the momentum and 
energy conservation, the relevant coupling term in H for 
this process will be: 



ftgk, q (b q al +q a k + h.c.) 



(7) 



Solving the amplitude equations. — As long as Brillouin 
cooling does not change the optical mode populations sig- 
nificantly, we can describe the interaction using coupled 
amplitude equations for the light field and the vibrational 
mode. We can write a Hamiltonian for the laser-driven 
mode cfi, the anti-Stokes mode d 2l and the phonon mode 
b of the form (written in a frame rotating at the laser 
frequency) : 
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This leads to the following classical equations for (a\) 
a\ etc.: 



di = [iA - «i/2] ai + K 1 /2a™ ax - iga 2 p* , 
d 2 = [i (A — Suj) — n 2 /2] a 2 — igot\P , 
p = [-in - T/2] p + £(t) - iga 2 a{ , 



(9) 
(10) 

(11) 



with (f *(*)£(*')) = S{t - t'). Here, af ax is the ampli- 
tude of the laser-driven mode at resonance (A = 0) in 
the absence of optomechanical coupling (g = 0). Su is 
the frequency difference between the optical modes. We 
solve this system analytically in the simplified case where 
the laser-drive mode is on resonance, and the frequency 



difference between modes is chosen to be equal to the 
mechanical mode frequency ft. Additionally, we assume 
a non-depleted pump and linearize the solution around 
(ol\,ol2,P) — (a7 iax , 0, 0). Taking these simplifications 
into account, the equations to be solved reduce to: 



d 2 = [-iScu - k 2 /2] a 2 - iga™ ax p 



P=[ift- T/2] P + vWf m - iga 2 a? 



(12) 
(13) 



We solve the system of equations in Fourier space to ob- 
tain an expression for (\p(t)\ 2 ) = n and (|o^2 (^) | 2 ) = n a Si 
the average phonon number and anti-Stokes photon oc- 
cupation respectively: 
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Quantum noise approach. — Alternatively, we can em- 
ploy the quantum noise approach to derive the cooling 
and heating rates for the mechanical mode q subject to 
a situation where the mode k is assumed to be laser- 
driven and photons are scattered into the other mode 
k + q. In contrast to the amplitude equations discussed 
above, the rate equation approach will only work for 
k » T opt + r&fj but unlike the amplitude approach it 
is not limited to situations where the optical mode pop- 
ulations remain essentially unchanged by the Brillouin 
processes. The idea will be to replace the driven photon 
mode operator by a c-number and to view the result- 
ing Hamiltonian as composed of a fluctuating quantum 
noise term coupling to the mechanical mode. This quan- 
tum noise term essentially arises from the interference 
between the driven mode and the vacuum fluctuations 
in the second mode. Consider the mode k to be driven 
. Upon substitution, the coupling term 



as a k = ae 
becomes: 



-iLJL,t 



hg k , q (b q a k+q ae lUht + h.c.) . 



(19) 



where b q couples to fluctuating quantum noise variable, 
F = hg kyq dl +q ae~ luJL \ in the form H int 
The transition rate for phonon annihilation is: 
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Likewise, the transition rate for phonon creation is: 

-Q 



(f t F) = J dte- iu>t (p1(t)F(0)) = . 



(22) 
(23) 



Taking into account both transition rates, we are left with 
a cooling rate that is set by the balance between up- and 
downward transitions. It can be written in the form: 



opt 



T pt — 9k+q 



(K/2)2 + (co k+q -LO L - ft) 2 ' 



(24) 
(25) 



where T opt quantifies the rate of Brillouin scattering. 
Here we have split off the dependence on the photon 
number n p hot circulating inside the lower optical mode. 
Below, we will display the slightly generalized expressions 
for the average up- and down-transition rates (deduced 
from Eqs. (20) and ([22]) ) , for the case of arbitrary photon 



numbers in both optical modes. 

Rate equation approach. — We can now proceed to 
solve a system of rate equations to determine the average 
phonon number in our system. We note that these rate 
equations do not take into account non-resonant scat- 
tering processes (scattering into the tails of the optical 
density of states, suppressed by a factor (n/Q) 2 ), and 
would also cease to be valid in a strong coupling regime 
(where T opt > Having derived the form of T opt: let 

us now consider a two level optical system with photon 
decay rates K\ and k,2. The transition rates between the 
two levels will be written in terms of the optomechanical 
coupling rate r , the average photon occupations of the 
two optical states n p , n a s, and the average phonon occu- 
pation n. The up transition corresponds to cooling and 
the down transition to heating, 



r t = T n p (n aS + l)n , 

r i = r (n p + l)n aS (n + l) . 



(26) 
(27) 



We can write the rate equations for the two levels based 
on the cavity decay rates and the heating and cooling 
transition rates, and the steady state photon number due 
to the driving laser , 



n a S = -n a S^2 ~ Y\ + I> , 

n v = ( n p ~ n P )Ki + r | - I> , 
n = (nth ~ n)Y M 



r+ - r t . 



(28) 
(29) 
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Solving for the steady state solution by setting the time 
derivatives equal to zero yields a set of equations which 
relate the photon occupations to the pump laser and 
acoustical phonon occupation: 

n p = rip - (r/«i) , (31) 
fi aS = (I> 2 ) , (32) 

r = r t - r+ = (n th - n)v M . (33) 
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FIG. 2. Average phonon number as a function of input power. 
The average phonon number starts from a thermal occupation 
of about 10 5 and cools by a ratio of about 260. Inset: Intra- 
cavity anti-Stokes dependence on input power. 



We see that Y = r[n,n p (n),n a s(n)] is a nonlinear func- 
tion of n. Thus we must solve the following relation, 



T(n) = (nth ~ n)T M • 



(34) 



where n will be a function of the dimensionless parame- 
ters n th ,n p , ^r, 

Discussion. — In the simple limit where the optical re- 
laxtion rates ^i, K2 are fast compared to the mechanical 
relaxation rate and optomechanical coupling rate r , Tm, 
we find only the cooling rate remains. In this case, the 
system simplifies to the following relations: 
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This is the usual classical cooling result where the cool- 
ing is unlimited with input power. We can also solve the 
full set of equations numerically, and obtain the cooling 
rate as a function of input power. For the system, we 
take for example a 100/im diameter Si02 sphere [2j [TO] 
such as the ones used in the context of Brillouin scatter- 
ing previously [8j |9]. The acoustical mode is taken to be 
a 50 MHz surface type mode like what was theoretically 
suggested [TT], numerically calculated [12 , and experi- 
mentally observed [9] to be excited via forward Brillouin 
scattering [13]. The sphere is taken to have two optical 
modes of high quality factor Q = 4xl0 8 . The optical 
modes are exactly separated by the mechanical oscilla- 
tion frequency. Our pump is a telecom compatible source 
(A = 1.55 jam). The relevant area is taken to be where the 
modes overlap, which is proportional to A 2 (the optical 
mode area) [14] as the optical mode is much smaller than 
the acoustical one [T2] . 

Figure [2] plots our theoretical prediction for the average 
phonon number in such a system as a function of input 
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FIG. 3. Average phonon number as a function of both in- 
put power and anti-Stokes quality factor. The dashed lines 
indicate the regime where our approximations hold. 



power. Significant cooling would begin at pump input 
powers of a few nW and saturates after a few mW . Start- 
ing from a room temperature phonon number of about 
10 5 (as is typical for such systems), one could cool by a 
ratio of about 260 in this example. As seen in the inset 
of Fig. [2j the amount of power circulating in the anti- 
Stokes mode clamps as the cooling process begins. As 
input power is increased, the final average phonon num- 
ber asymptotically converges to a lower limit: 



final 



n t h 



K>2 



(39) 



Eq. (39) indicates that if the quality factor of the anti- 



Stokes resonance were deliberately lowered compared 
to the pump resonance, higher cooling ratios could be 
acheived: Brillouin cooling requires an efficient way to 
get rid of the anti-Stokes photons. We therefore deter- 
mine the average phonon value as a function of both 
pump power and anti-Stokes mode quality factor, while 
all other parameters are left unchanged. Though many 
experimental challenges need to be overcome, as seen in 
Fig. |3j for diminished anti-Stokes quality factors, cooling 
ratios above 10 4 , can be acheived. 

Conclusion. — We describe here a triply resonant struc- 
ture for Brillouin cooling and develop the theory de- 
scribing this system. Unlike before [T5H2T] . for Brillouin 
cooling, the Doppler reflector is a monotonically travel- 
ling acoustic wave. Additionally, different from fluores- 
cent anti-Stokes cooling [22] that cools the whole thermal 
bath, Brillouin cooling evacuates heat from one selected 
vibrational mode, making it attractive for experiments 
in which this particular natural frequency is addressed 
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